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Radiative processes lead to important corrections to (e, e'p) experiments. While radiative correc- 
tions can be calculated exactly in QED and to a good accuracy also including hadronic corrections, 
these corrections cannot be included into data analyses to arbitrary orders exactly. Nevertheless 
consideration of multi-photon bremsstrahlung above the low-energy cut-off is important for many 
(e, e'p) experiments. To date, higher-order bremsstrahlung effects concerning electron scattering 
experiments have been implemented approximately by employing the soft-photon approximation 
(SPA). In this paper we propose a novel approach to multi-photon emission which partially removes 
the SPA from (e, e'p) experiments. In this combined approach one hard photon is treated exactly; 
and additional (softer) bremsstrahlung photons are taken into account resorting to the soft-photon 
approximation. This partial removal of the soft-photon approximation is shown to be relevant for 
the missing-energy distribution for several kinematic settings at mami and tjnaf energies. 



I. INTRODUCTION 

Coincidence electron scattering experiments are an 
important tool for probing both nuclear structure and 
the structure of the nucleons. A{e,e'p){A — 1) exper- 
iments {A being the atomic number) are e.g. used to 
study the independent particle shell model and spectral 
functions and, in connection with that, occupation 
numbers and correlations. And the (e, e'p) reaction has 
also been employed to study the structure of, e.g., the 
proton electric form factor in the so called Rosenbluth 
technique, apparently differing from polarization transfer 
measurements 

Results from electron scattering experiments are sub- 
ject to radiative corrections, i.e. QED amplitudes going 
beyond the leading-order Born term (see fig. [T|) . These 
corrections are of relative order a, but since they come 
with large logarithmic corrections they can contribute 
significantly to the cross section. Vertex correction, 
vacuum polarization, self-energy diagrams, and the 
two-photon exchange (TPE) are referred to as internal 
radiative corrections. And the four bremsstrahlung 
diagrams constitute the external radiative corrections 
and are the main focus of this paper. By introducing 
a small parameter associated with the photon energy 
resolution of the detectors, Eq, one can split up the 
cross section into a 'non-radiative part' including vertex 
corrections, vacuum polarization, self-energy contribu- 
tions, TPE, and the emission of soft bremsstrahlung 
photons with energies below Eq; and into a 'radiative 
part', accounting for the emission of bremsstrahlung 
photons with energies above the low-energy cut-off Eq 
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The individual contributions from the internal and 
external radiative correction diagrams are divergent. 
It has first been shown by Schwinger ^ that by in- 
troducing the low-energy cut-off these divergences cancel. 

The emission of bremsstrahlung alters the particle 
momenta and energies seen by the detectors and has 
to be corrected for in data analyses. Mo and Tsai first 
discussed this feature of electron scattering experiments 
0, 01 , considering single-photon bremsstrahlung exactly, 
aside from an approximation in the calculation of 
the TPE contribution. Multi-photon emission is only 
included for soft photons with energies smaller than Eq. 
A comprehensive review on radiative corrections can be 
found in ref. @. 

In order to obtain the desired experimental accu- 
racy radiative corrections cannot be limited to the 
second order amplitudes. One has to take into account 
higher-order bremsstrahlung processes (multi-photon 
emission) above the low-energy cut-off [1, 0, |§] which is 
straight-forward in soft-photon approximation (SPA). 
In SPA, multi-photon emission translates to a simple 
exponentiation of the bremsstrahlung contribution 
because of a factorisation of the bremsstrahlung cross 
section d, 0, Hi- A bremsstrahlung photon of energy 
uj'^ has no backlash on the electron-proton scattering in 
the limit where — > 0. In this limit bremsstrahlung 
amplitudes factorise into the first-order Born amplitude 
times the amplitude for emitting a soft photon. The 
factorisation also holds for multi-photon emission, if 
each emitted photon has vanishing energy. Summing 
over all orders of soft-photon emission then leads to 
exponentiation which also gives the cross section the 
correct asymptotic behaviour as the parameter Eq 
[E S Hi Si- Strictly speaking the SPA is only valid in 
the limit where all bremsstrahlung photon energies go 
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Feynman diagrams beyond the leading order. The 
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FIG. 1 

leading Feynman diagram A^B^m together with vertex cor 
rections, vacuum polarization and bremsstrahlung are called 
the Born approximation. 



to zero. In data analyses the approximation is applied 
however to emission of bremsstrahlung photons with 
finite energies [fOj. The question arises up to which 
upper bound for the maximum bremsstrahlung photon 
energy the SPA is a good approximation. 

While single-photon corrections to electron scattering 
experiments have been calculated exactly, including 
hadronic loops to a good accuracy 0, 0, Q and including 
proton structure [llj, multi-photon data analyses arc 
usually performed in SPA. It allows both for straight- 
forward inclusion of higher-order bremsstrahlung and 
for straight-forward Monte Carlo generation of the 
bremsstrahlung photon angular distribution 1Q\. The 
purpose of this paper is to improve the multi-photon 
radiative correction treatment for (e, e'p) experiments 
at MAMI and tjnaf energies by partially removing the 
SPA from multi-photon radiative corrections. 

As the dominant contribution to bremsstrahlung in 
(e, e'p) experiments comes from electron bremsstrahlung, 
we will initially omit bremsstrahlung originating from 
the proton. In sect. IIVI we will then include proton 
bremsstrahlung with recourse to SPA, showing that 
proton bremsstrahlung is only a minor correction to the 
full radiative corrections. 

The two TPE diagrams (see fig. [1]) go beyond the Born 



approximation and have received a lot of attention re- 
cently in H(e,e'») experiments studying the proton elec- 
tric form factor [l|, [12, [3] ■ But they are non-relevant for 
A(e,e'p) nuclear structure experiments where A> 1. 



II. MULTI-PHOTON SPA CROSS SECTION 

The emission of a single photon by, e.g., the scattered 
electron is 



(fc' -|- Lu)^ — m? 



x-v{p')rf,{q'')v{p), 



q 



(1) 



where k and k' are the four-momenta of the incident and 
the scattered electron; and p and p' are the four- momenta 
of the incident and final proton. The electron and pro- 
ton spinors are denoted by u and v, respectively. The 
momentum transfer to the proton is defined by 



k — k' — LU . 



(2) 



Here uj — ui^ (1,1, ft j) denotes the four-momentum of the 
bremsstrahlung photon and e from eq. ([T]) is its helicity. 
The vertex 
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contains the electromagnetic structure of the target nu- 
cleus. M and m are the proton and electron masses, 
respectively. In the limit where the photon energy uj^ 
vanishes the SPA can be applied and the single-photon 
emission amplitude llJ simplifies yielding [3,13] 
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Here the first-order Born amplitude is given by 
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The factorization of the SPA amplitude ([4]) also applies to 
the other three bremsstrahlung diagrams (see A^broms in 
Fig. [T]). Summing and averaging over the spins and per- 
forming the QED traces, we obtain the factorized single- 
photon cross section as 
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where 



is the elastic first-order electron-proton scat- 
tering Born cross section and 
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is the SPA angular distributions of the bremsstrahlung 
photon [8j]. 

This distribution ([7|) exhibits peaks as a function of 
the photon angles Q-y in the directions of the incident 
electron, the scattered electron, and the recoiling proton. 
Observing that most of the bremsstrahlung is emitted 
along the e-direction and the e'-direction, Schiff intro- 
duced the so-called peaking app roximation for inclusive 
(e, e') experiments in 1952 |16j| which was extended to 
exclusive {e,e'p) experimentsfsl. [TtI. [lH . approximating 
the angular distribution ([7]) by 
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and can be found e.g. in ref. [8|. The functions 
Ac, Ac', and Ap' can be obtained by integrating the 
e-, e'-, and p'-contribution to the angular distribu- 
tion ([7]). The peaking approximation can be removed 
from (e, e'p) data analyses, as has been shown in ref. p^ . 

Now we consider the SPA cross section for multi- 
photon bremsstrahlung with total photon energy below 
the low-energy cut-off Eg together with the emission of 
n photons with energies above the cut-off, which is given 
as M, 
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As mentioned above the exponential function accounts 
for multi-photon emission below the cut-off energy Eq to 
all orders. (^soft(£'o) (Shard contain the bremsstrah- 
lung and the internal radiative corrections, respectively. 
Their definitions can be found in ref. [1]. Necessary 
integration techniqes are given, e.g., in ref. [l9l |. 

In order to obtain an event generator we can sum the 
cross section ([9]) over all n photons and we integrate over 
the photon energies up to an upper limit E'tot , which is a 
parameter associated with the cuts applied to the data. 
We obtain the photon-integrated cross section (loj 
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where Xa is an " acceptance function" of the electron and 
photon kinematic variables (for those photons with en- 
ergies larger than Eq), Xa = Xa(^c, , . . . , uj", f2„), 
which is the probability for an event to be seen in [xa] 
[13] ■ In (fTO|) Etot is chosen such that all bremsstrahlung 
photons with non-vanishing Xa included in the inte- 
gration, and A is the integral of A{n~^) over the angular 
distribution, see Eg. ([7|), 
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Cross section ([T0| is expressed in terms of a probability 
density function which can be used to generate events 
with certain photon energies, multiplicities, and angular 
distributions using Monte Carlo event generators ^10| . 
Cross section (fTO|) is independent of the cut-off Eq in the 
limit where Eq becomes small (see Ref. |8|). 

The Monte Carlo routine we introduce in this pa- 
per generates multiple photons according to cross 
section (fTO|). The photon multiplicities n follow a 
Poisson distribution and the bremsstrahlun g en ergies 
are essentially distributed according to l/w^ [lO|. The 
bremsstrahlung photon angles are generated accord- 
ing to the elastic angular distribution Aciiill^) from 
eq. (O, using the elastic (e,e'p) kinematics for the 
electron and proton momenta, by a rejection algorithm. 
The index i indicates that each individual photon from 
a multi-photon event follows the angular distribution ([7]). 

The SPA simplifies multi-photon bremsstrahlung cal- 
culations considerably by the factorisation of the cross 
section, (e, e'p) data analysis procedures additionally 
adjust the kinematic settings of the {e,e'p) reaction. 
And the form factors are evaluated at the adjusted value 
of d, lisl . following a suggestion by Boric and 
Drechsel [20, UM- We refer to this extended version 
of the SPA as the 'modified SPA', or mSPA, which 
accounts for the bremsstrahlung photons' backlash on 
the electron-proton scattering process approximately via 
a modification of the kinematics. 



III. MULTI-PHOTON BREMSSTRAHLUNG 
EVENT GENERATION BEYOND SPA 

One way of extending the SPA for single-photon brems- 
strahlung is to calculate the cross section not in elastic 
kinematics but in the real 17-kinematic settings shown 
in fig. O Given the photon's four-momentum and the 
beam energy and fixing e.g. the electron scattering an- 
gles, one calculates the remaining kinematic quantities. 
The cross section is then calculated using the new 17- 
kinematic variables, additionally inserting the modified 
value of into the form factors. 
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FIG. 2; The I7 matrix element squared. 



Data analyses cannot include an exact multi-photon 
QED treatment, since the exact multi-photon amplitudes 
cannot be included into the calculations to arbitrarily 
high orders. Therefore we suggest a novel combined 
approach: we generate multi-photon bremsstrahlung 
according to the SPA distribution. We choose one 
hard photon from each scattering event which we treat 
exactly by calculating the exact QED I7 matrix element 
(see fig. H]) . And we treat the remaining photons as soft 
photons (with energies still larger than the low-energy 
cut-off £'0), employing the SPA. 

There are several approaches to choosing the hard pho- 
ton. We compare four different ones here: 

I. The photon with the largest energy tj" is chosen as 
the hard photon. This method addresses the ques- 
tion of the validity of the SPA. In this approach the 
choice for the hard photon depends on the choice 
of the reference frame, in contrast to the other ap- 
proaches described here. 



2. Calculating the momentum transfer squared, 
for each photon of the multi-photon event we 
choose the photon yielding a value of qf. 



which 



deviates most from the elastic value. 
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In this 



method form factors are evaluated using this qf^. 

3. Calculating the momentum transfer in multi- 
photon kinematics, q^^^, the photon is chosen which 
leads to a value of qf^ which is closest to g^ot . This 
method and method 2 focus on both the kinematic 
aspects of bremsstrahlung emission and on the ap- 
propriate calculation of the form factors. 

4. Choose the 'hard' photon randomly. This approach 
was tested in order to demonstrate the influence of 
the hard photons on the result. 

We now go beyond the SPA by re-weighting our SPA 
Monte Carlo generator described above, given by the 
ratio of the exact matrix element and the SPA matrix 
element. 

In order to simplify the type face, we absorb the photon 
energy dependence into the angular distribution A{ill^), 
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In SPA the 17-bremsstrahlung matrix element squared 
for the hard photon (see fig. [5]) then reads 
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Considering single-photon bremsstrahlung first, we as- 
sign to each Monte Carlo bremsstrahlung event a weight 
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This weight divides out the approximate SPA matrix el- 
ement and replaces it by the exact matrix element for 
single-photon radiation. In contrast to the 'exact weight' 
the mSPA weight measures the influence of the 
mSPA, 
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As the exact weight ([T4| this factor re- weights the event 
but now for the case of the SPA matrix elements in 
modified kinematics. So the first-order Born amplitude 
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as well as the photon distribution A„ispa(^) in 



the numerator are evaluated in 17-kinematics (or mSPA 
kinematics which is the same in this case). 

This approach can be extended to multi-photon brems- 
strahlung. Assuming that the n bremsstrahlung photons 
have been ordered such that the nth photon is the hard 
photon, we define the exact weight as 
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where AmSPA is the angular photon distribution in mod- 
ified SPA; and we define the modified weight as 
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(17) 

The Monte Carlo routine introduced here generates 
multi-photon events in SPA. Each scattering event is as- 
signed two weights and p?)) using the four methods 
for selecting the hard photon, described at the beginning 
of this section. In order to check the combined approach 
against data (where available) and against SPA simula- 
tions, we embedded the combined Monte Carlo routine 
into SIMC, a Monte Carlo simulation for electron scatter- 
ing experiments from tjnaf [i3|. This data analysis 
code uses a version of the mSPA which - in contrast 
to our mSPA calculations shown in figs. [3] to [5] - is not 
able to choose a hard photon from a given multi-photon 
bremsstrahlung event. On top of that it makes use of 
the peaking approximation fluL \it\ . Thus we were able 
to compare the error associated with the use of the SPA 
with other sources of errors. 
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FIG. 3: Missing-energy distribution for multi-photon 
bremsstrahlung. The solid curve (1) (green) was obtained us- 
ing the modified SPA weight ffTT} . The dotted curve (2) (red) 
represents the combined approach using weight ([16} . The in- 
set graph shows the deviation between the two curves in per- 
cent normalized to the combined result. At Em = 100 MeV 
the modified SPA calculation overestimates the radiative tail 
by 3.6%, at = 500 MeV the deviation is 25%. The mo- 
mentum transfer is — 0.6 GeV^. 
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FIG. 4: Missing energy distribution for multi-photon 
bremsstrahlung. The attribution of the curves is as above in 
fig. El At Sni = 100 MeV the modified SPA calculation over- 
estimates the radiative tail by 3.9%, at Em ~ 500 MeV the 
deviation is 29%. The momentum transfer is — 2.0 GeV^. 



IV. RESULTS 

To test the combined calculation we simulated the 
missing-energy distributions and the photon angular 
distributions at several kinematic settings. Ten million 
events per run were generated in order to get good 
statistics for the different yields. The results turned out 
to be indistinguishable under approaches 1, 2, and 3 for 
choosing the hard photon (see previous section). Only 
the random choice of the hard photon (approach 4) 
deviated from the other methods. The results presented 
in this section are therefore based on approach 1. 
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FIG. 5: Missing-energy distribution for multi-photon 
bremsstrahlung. The attributions of the curves is as above in 
fig. El At Bni = 100 MeV the modified SPA calculation over- 
estimates the radiative tail by 3.5%, at Em = 500 MeV the 
deviation is 23%. The momentum transfer is — 4.0 GeV^. 
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FIG. 6: Missing-energy distribution for multi-photon 
bremsstrahlung. The attribution of the curves is as above in 
fig. E] At Em = 100 MeV the modified SPA calculation over- 
estimates the radiative tail by 3.0%, at Em = 500 MeV the 
deviation is 17%. The momentum transfer is — 8.0 GeV^. 



Figures. [3] to [S] show multi-photon missing-energy 
distributions, once calculated in mSPA P?|) . once cal- 
culated using the combined approach (fT6|) . We see that 
the mSPA calculations overestimate the radiative tails. 
While the deviations between the two calculations are 
of the order of a few percent for missing energies below 
= 100 MeV, the deviations become considerably 
larger towards the far ends of the radiative tails. 

Since we want to evaluate the importance of the SPA 
in multi-photon radiative corrections it is important to 
compare the effect of the improved radiative corrections 
to other sources of errors, present in (e, e'p) experiments. 
On top of radiative corrections, (e, e'p) experiments are 
e.g. corrected for finite detector resolution, for particle 
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FIG. 7: Missing-energy distribution for multi-photon 
bremsstrahlung simulated with SIMC. The solid curve (1) 
(black) represents the standard SiMC modified SPA radiative 
corrections. The dotted curve (2) (red) shows the Em distri- 
bution obtained by inserting our combined radiative correc- 
tion approach into SiMC. The latter one has more strength 
in the radiative tail. The total yield differs by 4.6%. The 
momentum transfer is = 0.6 GeV^, as in fig. [S] 
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FIG. 8: Missing-energy distribution for multi-photon 
bremsstrahlung simulated with SIMC. The attribution of the 
curves is the same as in fig. [3 The total yield differs by 2.3%. 
The momentum transfer is — 2.0 GeV^, as in fig. U 

decays, and for multiple scattering pTl-ITsj. 

The simulations shown in figs. [7| through [TO] take 
into account radiative corrections, detector resolution, 
particle decays, and multiple scattering. They indicate 
that the combined approach has an impact on the 
missing-energy distribution which is not diluted by 
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FIG. 9: Missing-energy distribution for multi-photon 
bremsstrahlung simulated with SIMC. The attribution of the 
curves is the same as in fig.[7| The total yield differs by 1.4%. 
The momentum transfer is = 4.0 GeV^, as in fig. [5| 
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FIG. 10: Missing-energy distribution for multi-photon 
bremsstrahlung simulated with SIMC. The attribution of the 
curves is the same as in fig. [7] The total yield differs by 1.4%. 
The momentum transfer is = 8.0 GeV^, as in fig. |6| 



the other corrections. The radiative tails calculated 
following the combined approach are stronger than 
the radiative tails obtained with the standard SIMC 
radiative correction procedure for all kinematic settings 
considered. 

As a second observable we considered the angular 
distribution of the bremsstrahlung photons. As for the 
missing-energy distributions we first looked at radiative 
corrections only. Then we additionally took into account 
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FIG. 11: Angular bremsstrahlung photon distribution. The 
solid curve (1) (green) represents the modified SPA, the dot- 
ted curve (2) (red) shows the combined approach. The SPA 
distribution deviates from the combined calculation especially 
in the vicinity of the incident electron, around 0°. The kine- 
matic settings are as in fig. O 
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FIG. 12: Angular bremsstrahlung photon distribution. The 
attribution of the curves is as in fig. 1111 The SPA distribu- 
tion deviates from the combined calculation especially in the 
vicinity of the incident electron, around 0°. The kinematic 
settings are as in fig. |31 



FIG. 13: Angular bremsstrahlung photon distribution. The 
attribution of the curves is as in fig. 1111 The SPA distribu- 
tion deviates from the combined calculation especially in the 
vicinity of the incident electron, around 0°. The kinematic 
settings are as in fig. [5] 
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FIG. 14: Angular bremsstrahlung photon distribution. The 
attribution of the curves is as in fig. 1111 The SPA distribu- 
tion deviates from the combined calculation especially in the 
vicinity of the incident electron, around 0°. The kinematic 
settings are as in fig. |6] 



other corrections. The photon angular distributions are 
shown in figs. [TT] to [TH For aU Icinematic settings the 
mSPA calculations overestimate the angular distribution 
in the vicinity of the incident electron, that is small 
angles. As for the missing energy the question is whether 
this deviation can be seen in the data, or whether other 
sources of errors dominate the photon angular distribu- 
tion. In order to determine the impact of the combined 
approach on the photon angular distribution we resort 
to SIMC, as in the case of the missing-energy distribution. 

Using standard SiMC for comparison implies that the 
peaking approximation is employed [iTj . It is clear that 
the photon angular distribution, generated in SiMC's 
standard radiative correction approach differs from the 



combined approach even more than the mSPA from 
the combined calculation shown in figs. [TT] to [T3] (see 
also ref. (lo|'). The peaking approximation and thus the 
standard SiMC analysis code is known not to describe 
the experimental angular distribution accurately in 
between the two radiation peaks coming from the 
incident and the scattered electron fTo'|. Nevertheless, 
looking at the photon angular distributions both with 
the standard SIMC code and with the modified version 
of SIMC, containing our combined simulation, gives us 
the opportunity to rate the impact of our approach 
compared to the standard SIMC radiative correction 
procedure. Where available, we compared the different 
approaches to data. 

Figs. [15] to [TH] show the photon angular distributions 
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FIG. 15: Angular bremsstrahlung photon distribution for 
H(e,e'p) scattering generated with SIMC. The sohd curve (1) 
(black) represents the standard SiMC modified SPA calcula- 
tion. The dotted curve (2) (red) shows the angular distri- 
bution simulated with the combined approach. The largest 
difference between the standard SiMC treatment and the data 
occurs in the middle between the e and the e' directions and 
in the height of the two peaks. This is due to the peaking 
approximation which is used by the standard SIMC code; it 
overestimates the peaks. The kinematics are the same as in 
figs. Eland [Til 



as generated by SIMC. The largest deviations between 
the standard SIMC radiative correction and the combined 
approach appear in the middle between the peaks due to 
e and eJ bremsstrahlmig. This is in contrast to figs. [11] 
to [T4| where deviations occur at small angles. This 
indicates that the small deviations between SPA and 
combined approach seen at small angles (see figs. [11] to 
[Hj) are washed out by other corrections. The photon 
angular distribution is not effected by the SPA and by 
the combined approach. The dominant approximation 
is the peaking approximation here. 

Fig.[T9]shows experimental data for the photon angular 
distribution reconstructed from the H(e, e'p) kinematics 
shown in figs.[3][Tll and [T5l (tjnaf experiment E97-006, 
[2^ [23II) and verifies this conclusion. This measurement 
at = 2.0 GeV'^ shows that both the mSPA calculation 
and the combined approach describe the data well. 

In order to additionally include the effect of proton 
bremsstrahlung in the combined approach, we inserted it, 
assuming the mSPA is a good approximation for proton 
bremsstrahlung, the overall effect of proton bremsstrah- 
lung being small anyway. The first-order electron brems- 
strahlung is still calculated exactly. This was achieved 
by modifying the weight ([TO]) into 
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FIG. 16: Angular bremsstrahlung photon distribution gener- 
ated with SIMC. The attribution of the curves is the same 
as in fig. 1151 The comparison with experimental data can be 
found in fig. [I9l The kinematics are the same as in figs. |4]and 
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FIG. 17: Angular bremsstrahlung photon distribution gener- 
ated with SIMC. The attribution of the curves is the same as 
in fig. [15] but there was no data available for these kinematic 
settings. They are the same as in figs. [5l and [131 
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where Arcst(whard) is the angular distribution ([7]) 
without the electron-electron terms. (In addition ([7]) 
has been divided by We generated the photon 

angular distribution using this approach, comparing it 
to the combined calculation without the proton (|16p . 

As can be seen in fig. [^D] including mSPA proton 
bremsstrahlung does not change the photon angular 
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FIG. 18: Angular bremsstrahlung photon distribution gener- 
ated with SIMC. The attribution of the curves is the same as 
in fig. [15] but there was no data available for these kinematic 
settings. They are the same as in figs. 151 and [T^ 



distribution considerably. Just the electron radiation 
peaks are slightly overestimated when neglecting the 
proton bremsstrahlung. The missing energy is not 
changed at all by including proton bremsstrahlung (see 
fig. [21]). We conclude that the proton bremsstrahlung is 
not relevant and that mSPA is a good approximation for 
proton bremsstrahlung. 

We observe that in the case considered here the com- 
bined approach reproduces the experimental photon an- 
gular distribution much better than the standard SiMC 
simulation using the peaking approximation. The com- 
bined approach presented in this manuscript and the full 
angular approach from ref. [lo| are much closer to each 
other and to the data than the standard SiMC simula- 
tions, indicating, again, that removal of the peaking ap- 
proximation is more important for the photon angular 
distribution than the partial removal of the SPA. 



V. DISCUSSION AND CONCLUSION 

We have shown that it is feasible to improve (e, e'p) 
radiative corrections by partially removing the SPA from 
multi-photon bremsstrahlung processes. Our results are 
invariant under three different approches to choosing the 
hard photon whereas a random choice of the hard photon 
is not feasible. Treating one photon as a hard photon 
improves the kinematic treatment of the (e, e'p) reaction. 

We compare the combined approach to multi-photon 
mSPA calculations, looking at the missing-energy distri- 
bution. The mSPA overestimates the radiative tail for 
different kinematic settings. In order to check whether 
the combined treatment has an influence when consider- 
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FIG. 19; Angular bremsstrahlung photon distribution for 
H(e,e'p) scattering generated with SIMC and compared to 
data. The bump on the right hand side at the proton di- 
rection is an artefact due to 'punch through effects'. The 
dash-dotted curve (1) (green) and the (2) (red) curve show 
standard SiMC photon distributions, both using different ver- 
sions of the peaking approximation. The dotted curve (3) 
(black) represents the data. The solid line (4) (magenta) is 
the new combined approach implemented into SIMC. Addi- 
tionally there is a blue line (5) from ref. [13] showing a full 
angular SPA simulation. The kinematics are the same as in 
figs. Ill and 1121 The combined approach (4) fits the data well. 
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FIG. 20: Photon angular distribution with (1) (red curve) 
and without (2) (black curve) SPA proton bremsstrahlung. 
The simulation neglecting proton bremsstrahlung overesti- 
mates the electron peaks slightly. In the rest of the photon 
angle domain the two curves coincide. 



ing additional experimental corrections we subsequently 
inserted our combined calculation into an existing data 
analysis code, which in its standard version is using 
the peaking approximation and a version of mSPA, in 
addition to other corrections [1, [l3|. We showed that 
that our combined approach had an impact on the 
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FIG. 21: Missing energy with (1) (red curve) and without (2) 
(black curve) SPA proton bremsstrahlung. The two curves 
are hardly distinguishable. 

missing-energy distribution and that it was not hidden 
by other sources of background and other effects like, 
e.g. detector resolution. The computational expense of 
the combined method was small, at most a factor of 2, 
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